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Abstract 

We estimate the mean first time, called the mean rotation time 
(MRT), for a planar random polymer to wind around a point. This 
polymer is modeled as a collection of n rods, each of them being pa- 
rameterized by a Brownian angle. We are led to study the sum of 
i.i.d. imaginary exponentials with one dimensional Brownian motions 
as arguments. We find that the free end of the polymer satisfies a 
novel stochastic equation with a nonlinear time function. Finally, we 
obtain an asymptotic formula for the MRT, whose leading order term 
depends on y/n and, interestingly, depends weakly on the mean ini- 
tial configuration. Our analytical results are confirmed by Brownian 
simulations. 

1 Introduction 

This paper focusses mainly on some properties of a planar polymer motion 
and in particular, on the mean time that a rotation is completed around a 
fixed point. This mean rotation time (MRT) provides a quantification for the 
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transition between a free two dimensional Brownian motion and a restricted 
motion. To study the distribution of this time, we shall use a simplified 
model of a polymer made of a collection of n two dimensional connected rods 
with Brownian random angles (Figure [1]). To estimate the mean rotation 
time, we shall fix one end of the polymer. We shall examine how the MRT 
depends on various parameters such as the diffusion constant, the number 
of rods or their common length. Using some approximations and numerical 
simulations, the mean time for the two polymer ends to meet was estimated 
in dimension three |WiF74t IPZS96j . 

Although the windings of a (planar) Brownian motion, that is the number 
of rotations around one point in dimension 2 or a line in dimension three and 
its asymptotic behavior have been studied quite extensively |Spi58[ IPiY86t 



ILeGY87] ILeG90[ IReY99] . little seems to be known about the mean time for 
a rotation to be completed for the first time. For an Ornstein-Uhlenbeck 
process, the first time that it hits the boundary of a given cone was recently 
estimated |VaklO] . 

The paper is organized as follows: in section 2, we present the polymer 
model. In section 3, we study a sum of i.i.d. exponentials of one dimensional 
Brownian motions. Interestingly, using a Central Limit Theorem, we obtain 
a new stochastic equation for the limit process. This equation describes the 
motion of the free polymer end. In section 4, we obtain our main result which 
is an asymptotic formula for E [r„] the MRT when the polymer is made of 
n rods of equal lengths Zq with the first end fixed at a distance L from the 
origin and the Brownian motion is characterized by its rotational diffusion 
constant D. We find that the MRT depends logarithmically on the mean 
initial configuration and for uIq » L and n > 3, the leading order term is 
given by: 

1. for a general initial configuration: 



9fc(0) 



21" I ^ Ti ^^'^ 



where (^^(0), 1 <k < n) are the initial angles of the polymer and Q 
9.56, 

2. for an initially stretched polymer: 



E [Tn] ^ ^ (ln(n) + Q) 
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3. for an average over uniformly distributed initial angles: 



E \tJ ^ — Q, 



where Q ^ 9.62. 

We confirm our analytical results with generic Brownian simulations with 
normalized parameters . Finally, in section 5 we discuss some related open 
questions. 

2 Stochastic modeling of a planar polymer 

Various models are available to study polymers: the Rouse model con- 
sists of a collection of beads connected by springs, while more sophisti- 
cated models account for bending, torsion and specific mechanical properties 
|Rou53trSSS80tlDoE94] . Here, we consider a very crude approximation where 
a planar polymer is modeled as a collection of n rigid rods, with equal fixed 
lengths Iq and we denote their extremities by {Xq, Xi, X2, . . . , X„) (Figured]) 
in a framework with origin 0. We shall fix one of the polymer ends Xq = {L, 0) 
(where L > 0) on the x-axis. The dynamics of the i-th rod is characterized 
by its angle 6i{t) with respect to the x-axis. The overall polymer dynam- 
ics is thus characterized by the angles {Oi{t), 02{t), . . . ,6n{t),t > 0). Due to 
the thermal collisions in the medium, each angle follows a Brownian motion. 

Thus, with ^^=^ denoting equality in law. 



{ei{t),i<n) =' V2D{Bi{t),i<n) that is: d6i{t) = V2D dBi{t), i < n, 



where D is the rotational diffusion constant and {Bi{t), . . . ,Bn{t),t > 0) is 
an n-dimensional Brownian motion (BM). The position of each rod can now 
be obtained as: 



X,{t) = L + /oe^'^W, X^it) = X,{t) + /oe^'^W, . . . , X„(t) = X„_i(t) + loe''"^'\ 

(1) 
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Figure 1: Schematic representation of a planar polymer winding 
around the origin, (a) A random configuration, (b) MRT when the n-th 
bead reaches 2tt. 



In particular, the moving end is given by: 



n 



n 



,iV2DBkit) 



k=l k=l 



(2) 



which can be written as: 



X„(t) = Rnit)e 



(3) 



Thus ifn{t) accounts for the rotation of the polymer with respect to the origin 
and Rn is the distance to the origin. 

In order to compute the MRT, we shall study a sum of exponentials 
of Brownian motions, a topic which often leads to surprising computations 
YorOl] . First, we scale the space and time variables as follows: 



Equation ([2]) becomes: 



n 




(5) 



k=l 
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where (^Bi{t), . . . , Bn{t),t > 0^ is an n-dimensional Brownian motion (BM), 
and for A; = 1, . . . , n, using the scahng property of Brownian motion, we have: 

Before describing our approach, we first discuss the mean initial configu- 
ration of the polymer. It is given by: 

where the initial angles ^a;(0) are such that the polymer has not already made 
a loop. After scaling, the mean initial configuration becomes: 



with 

4(0) = ^fe(0). (8) 

Prom now on, we use 9k instead of ^fe, B instead of B and t instead of i. 
Any segment in the interior of the polymer can hit the angle 27r around the 
origin, but we will not consider this as a winding event, although we could 
and in that case, the MRT would be different. Rather, we shall only consider 
that, given an initial configuration c„, the MRT is defined as: 

MRT=E{Tn\Cn} = E[Tn], (9) 

where 

T„ = inf{t >0,|<^„(t)| = 27r}. (10) 

Thus, an initial configuration is not winding when 

bn(0)|<27r. (11) 

Then, we can define the winding event using a one dimensional variable only. 
In general, winding is a rare event and we expect that the MRT will depend 
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crucially on the length of the polymer which will be quite long. Interestingly, 
the rotation is accomplished when the angle ifn{t) reaches 27t or — 27r, but 
the distance of the free end point to the origin is not fixed, leading to a 
one dimensional free parameter space. This undefined position is in favor of 
a winding time that is not too large when compared to any narrow escape 
problem where a Brownian particle has to find a small target in a confined 
domain |WaK93l IHoS04l ISSHE061 ISSH071 IBKH071 IPWPSK09] . 

In this study, we consider not only that the initial condition satisfies 
|V5n(0)| < 27r, but we impose that v?n(0) is located far enough from 27r, to 
avoid studying any boundary layer effect, which would lead to a different 
MRT law. Indeed, starting inside the boundary layer for a narrow escape 
type problem leads to specific escape laws |SSHE06] . Given a small £ > we 
shall consider the space of configurations such that |(y9„(0)| < 27i — e. We 
shall mainly focus on the stretched polymer, 

(^i(0),^2(0),...,^„(0)) = (0,0,...,0), (12) 

and thus c„ = n > (in this case, (pn{0) = 0). Finally, it is quite obvious 
that winding occurs only when the condition: 

nlo > L (13) 

is satisfied, which we assume all along. 

The outline of our method is: first we show that the sum X„(t) converges 
(eq. ([2])) and we obtain a Central Limit Theorem. Using Ito calculus, we 
study the sequence: 

-^Xn{t) = [Xn{t) - E (X„(t))] . (14) 



and prove that ^X„()f:), for n large, converges to a stochastic process which 
is a generalization of an Ornstein-Uhlenbeck process (GOUP), containing a 
time dependent deterministic drift c„e~*. This GOUP is driven by a mar- 
tingale [m^^ , t > j that we further characterize. Interestingly, the two 

cartesian coordinates of ^M^*-"'' , t > j converge to two independent Brown- 
ian motions with two different time scale functions. To obtain an asymptotic 
formula for the MRT, we show that in the long time asymptotics, where 
winding occurs, the GOUP can be approximated by a standard Ornstein- 
Uhlenbeck process (OUP). Using some properties of the GOUP |VaklO] . 
we finally derive the MRT for the polymer which is the mean time that 
|v?„(t)| = 27r. 
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3 Properties of the free polymer end Xn{t) us- 
ing a Central Limit Theorem 



In this section we study some properties of the free polymer end In 
particular, using a Central Limit Theorem, we show that the limit process 
satisfies a stochastic equation of a new type. To study the random part of 
Xn{t), we shall remove from it its first moment and we shall now consider 
the asymptotic behavior of the drift-less sequence: 



We start by computing the first moment E{Xn{t)). Because (see (|2])), 
[9i{t),t > 0) are assumed to be n independent identically distributed (iid) 
Brownian motions with variance 2D, after rescaling, we obtain: 




:X„(t) 



1 



[X„(t)-E(X„(t))]. 



(15) 



n 



fc=i 




(16) 



where 



Cn is defined by ([7]) and we have used that: 



E ^gi(Bfc(t)-Bft(0))j 

We study the sequence (fTSl) as follows: 



= e 2 



(17) 




(18) 



where F^it) 



giBfc(t) _ (^e'^^kit)y Applying Ito's formula to 




n 



(19) 
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with Zq^^ = 0, we obtain: 



= M,(")-i / (20) 



where: 



= ^/ V(^cos(5fc(s))-sin(5fc(s))) 
V ^ Jo 



k=l 

= -Sl^'^+tCtl (21) 

We shall now study the asymptotic limit of the martingales M^""^ as n — )■ 
oo and summarize our result in the following theorem; the convergence in 
law being considered there is associated with the topology of the uniform 
convergence on compact sets of the functions in C(l 



[)2\ 



Theorem 3.1 The sequence converges in law to a Brownian 

motion in dimension 2, with 2 deterministic time changes. More precisely, 

(St\Ct\t> 0) m (<.(ij^.^„..-,.,),7(jj;;,,„,.-.,),i > 0). (22) 

where {au,'ju,u > 0) are two independent Brownian motions. 

Remark 3.2 a) The convergence in law l[2^) is a new result and it shows 
that the sum of the complex exponentials of i.i.d. Brownian motions can be 
approximated by a two dimensional Brownian motion, with a different time 
scale for each coordinate. 

b) An extension of Theorem \3.1\ when (exp {iB{s)) , s > 0) a Brownian mo- 
tion on the unit circle is replaced by a BM on the unit sphere in M" is obtained 
in JHVYlll 

Proof of Theorem 13.11 See Appendix |X1 ■ 
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We conclude that the sequence converges in law: 

{Mt\t > 0) X +27A,i-e--v^ > 0). (23) 

The process {z'f'\t > 0), which we defined in (I20p . is a generalization of the 
classical Ornstein-Uhlenbeck process; it is driven by (M^*-"\t > 0) and, from 
(l20l). we obtain: 



ft 



Zf"^ = e-M etdMf). (24) 







Corollary 3.3 a) The sequence {z[^\t > 0) converges in law: 

^(n) ila^ ^(oo)^ (25) 



with: 



Z^(°^) = e-5 /" v/sinh(s) d6, + i e's /" v/cosh(s) (26) 
Jo Jo 



where ySt,5t,t > Oj are two independent 1- dimensional Brownian motions. 

b) Z[^^ ^-^^^ 4= ( + iN] , where N and N are two centered and reduced 
Gaussian variables. 

Proof of Corollary |33] 

a) From (EH), (ESD and (EH) we deduce: 



Zf^ = e--2 TeidMi") (27) 







n— >oo 



e 2 / 62 W ^ d5s + ie 2 / 62 W d5s = ZrK 

Jo ^ ^ , Jo ^ ^ ^ 

sinh(s) -y/ cosh(s) 

(28) 
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b) From fl28p . we change variables: u = t — s and we obtain: 







* _u 1 e-2(t-«) /•* _„ /l e-2(t-«) ^ 
e M/ 7:: d6u + i e 2 W - H d6^ 



u=t—s 

(law) Jr. V 2 2 V 2 



(29) 

e 2 c?,5^ + i —— / 6 2 c?,5^ ^ (30) 



-\/2 Jo "\/2 Jo 



where the two variables on the RHS of f l30l) are centered Gaussian with 
variance 1/2 and the convergence in for t — )• 00 may be proved by using 
the dominated convergence theorem. ■ 

Remark 3.4 In flt683^ and jOchiSS^ there is a more complete asymptotic 
study for z[^\ip), as n 00, for a large class of functions (p, hut we do not 
pursue this line of research here. 



On one hand, from the identities (ITSj) . (fT6l) . ( ITSl) and ( IT9l) . we obtain the 
following expansion: 

X„(t) = Xn{t)+E[Xn{t)] 

^Xn{t) = X„(t)-/'= v^zf) + c„e-^ (31) 

with n the number of rods/beads, zj:""^ a GOUP driven by M^"'^ which is 
given by (123|) . Cn = E {Y^'^=i e*^'='-°'') a constant depending on the mean initial 
configuration and / = ^ the rescaled distance of the fixed end from the origin 
{Iq is the fixed length of the rods). 

On the other hand, with Corollary 13.3( a). we deduce that: 

zt^ ^ ^X„(t)^^[x„(t)-E(x„(t))] 



n \/n 

1 



{law) 



e 2 

n— >-cxD 



Xn{t) - CnC 2 - / 

3 f v/sinh(s) dSs + 26-5 ! v/cosh(s) d~5s = Zl°^\ (32) 
Jo Jo 
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4 Asymptotic expression for the MRT 

We now study more precisely the different time scales of the two Brownian 
motions in (123|) and ( 132|) . We estimate Zf°°^ for t large, the regime for which 
the rotation will be accomplished. 

We introduce here the following notation: ^ denotes closeness in the L^— norm: 



Wf'^^ ^ Wf"^^ means that lim E 

t—^oo 



for two stochastic processes {Wt^^\t > 0) and > 0), the notation 

0. 



We shall show that, with (Mt = 5t + i5t,t > a 2-dimensional Brownian 
motion starting from 1: 

zj.°°^ = e~2 / y^sinh(s) + z e"2 / a/ cosh(s) d6s 
Jo Jo 



L2 t /"* 62 



e-2 / —dM,. (33) 
^0 V 2 

For this, it suffices to use the expression (l30l) and the following Proposition, 
which reinforces the ^ result in fl33D. 



Proposition 4.1 As t ^ oo, the Gaussian martingales 



(^j^ v/sinh(s) d5s- d6s,t> 0^ 



and 



/ A/cosh(s) — / —= d6s,t > 
Jo Jo v2 



converge a.s. and in L^. T/ie /zmzt variables are Gaussian with variances 
and —1 + 2-\/2 — 2as(l) ^ 0,033, where as{x) = argsinh(x) = log(x + 
a/1 + x^), a; G M, respectively. 

Thus, by multiplying both processes by e~^, we obtain ( !33|) . 

Proof of Proposition l4TT] The Gaussian martingale /q ^A/sinh(s) — 
has increasing process 



e-2« - 1 rfs 
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which converges as t — )■ oo. Hence, the hmit variable smh(s) — dSs 

is Gaussian, and its variance is given by (we change variables: u = e"^'' and 
B{a, b) denotes the Beta function with arguments a and ^): 



ds e"" 



-2s 



u 



[ du u-^/^ ((1 -u)- 2^1^ + 1) 
Jo 



1 



2B 



1 3 
2' 2 



7r-3 



To be rigorous, the integral jj^ du — u — l) ^, which is well defined for 

< a < 1, can be extended analytically for any complex a with Re(a) < 3. 
For the convergence of the second process, it suffices to replace sinh(s) by 

cosh(s) = y(1 + e~^*). The limit variable ^A/cosh(s) — d5s is 

also Gaussian, and, repeating the previous calculation, we easily compute its 
variance. ■ 



Asymptotic expression for the MRT 

For t large, we derive an asymptotic value for the MRT. First, from ( l24l) and 
(BTl): 



where, the sequence c„ is: 



(34) 



n 



(35) 



n is the number of rods/beads, and Cn = E (^^=1 e*^* *■"■') is a constant de- 
pending on the mean initial configuration. Thus, from (l3Tll . ( I33ll . ( IMll and 
also using the scaling property of Brownian motion: 



x„(t) ^'-^ y^y/") 



(36) 



*We recall that if {T{x),x > 0) denotes the Gamma function, then B{a,b) — ^-^{l+t) • 
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where 

r/"^ =e'i[cn + I eirfl,/2 ) . (37) 



Changing time and expression of the MRT 

To express the MRT, we now apply deterministic time changes. To make our 
writing simple, we denote Yt for Y^^\ and changing variables m = | in fl37|) . 
we obtain: 

Y2t = e-' (cn + [ e^M^) . (38) 



\ ^0 / 

Now, there is another BM {Mt,t > 0), starting from c^, such that: 

Y2t = e'' (b,,) , (39) 

where: 

at = I e^'^ds 



2Me - - i 

2 ' 



hence: 



a-^(t) = -ln(l + 2t). (40) 



2 

Applying Ito's formula to fl39|) . we obtain: 

We divide by and we obtain: 

dY2s , dl 

Thus: 



, ,dY2s\ , Id! 
im = im 



Y 



2s 
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which means that, if we denote: 



ef ^ Im( 







,t > 0, 



(41) 



the continuous winding process associated to a generic stochastic process Z, 
then: 

Thus, the first hitting times of the symmetric conic boundary of angle c 

Tl^""! =inf {t >0: l^fl =c}, (42) 

and 



Tr=mfU>0: 



Of 



(43) 



for an Ornstein-Uhlenbeck process Y with parameter A = 1, and for a Brow- 
nian motion B respectively, with relation (HOl) . satisfy: 



2Tri = ^ln(l + 2rfl). 



Finally, 



E 



2T\f\ 



and equivalent ly: 



E 



-E 



In 2 1 

+ -E 

2 2 



ln2 1 

4 4 



In ( 1 + 2Tf\ 



In ( Tf\ + I 



(44) 

(45) 
(46) 



In Tf\ + 



(47) 



Thus, by taking c = 27r, for n large, with (pn{t) denoting the total angle of 
X„(t), the mean time E [f„], where f„ = inf{t > 0, |<^n(^)| = 27r}, that X„(t) 
rotates around 0, is: 



\n2 + E 



m ( Tifi + 1 



(48) 
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By using the series expansion of log(l + x), we obtain informally: 



E 



In f Tf I + ^ 



E 



In rj^"i 



E 



In 1 + 



2T'' 



(49) 



-E 



-E 



.T 



+ 



(50) 



where this expansion should be understood as follows: we obtain an alter- 
nate series of negative moments of tJ^ ', which converges |VaYll| . For our 
purpose, we will truncate the series. We shall estimate and use only the first 

moment of 1/Tc . By numerical calculations, the other moments can be 
neglected. 

From the skew product representation |KaSh88i IReY99] there is another 
planar Brownian motion (/?„ + ijujU > 0), starting from logc„ + iO, such 
that: 



(/3« + hu) 



u=Ht=j' 



log 

and equivalent ly: 

Thus, with tJ^' = inf {t > : \-ft\ = c} and T^"^' = inf <{ t > : 
because 6 



log 



l3Ht; Ot = -fHf 



(51) 



(52) 



|9B[ — iH 



hence tJ^ ' = if,. 



n|7l' 



where: 



^ = inf{t : Ht> u} = / dsexp{2/3s) = A^, 



and for u = t]'^' , we obtain: 



rp\J 



(53) 



(54) 
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It is of interest for our purpose to look at the first negative moment of 
At which has the following integral representation for any t > ( jPufOO] . 
|D-MMY00j . p.49, Prop. 7, formula (15)): 



E 



' 1 " 









^"ye-^'*/2^oth(|y) dy. 



(55) 



With /3, = log(c„) + pf\ where (/3f , s > 0) is a one-dimensional Brownian 
motion starting from 0, we obtain: 

Tfl = H-\T}?\)^ / ds exp(2/3,) 



b(i), 



ds exp (2/3f 



where tJ^ ' 



= inf <^ t > : 



1) 



(56) 



c| is the first hitting time of the sym- 



metric conic boundary of angle c of a Brownian motion B*^^) starting from 
1+iO. Hence, from (1541155^ and by using the Laplace transform for the hitting 
time T^' |ReY99] (Chapter II, Prop. 3.7) or |PiY03](p.298): 



E 



e 2 



1 



cosh(?/c) 



(57) 



we get: 



E 



^|gB(l)| 



yE 



y 



cosh(?/c) 

For c = 27r, we obtain the numerical result: 



G(27r) = E 



coth (^^y^ dy 
coth (^^y^ dy = G{c). 



b{i)i 



27T 



0.167. 



Thus, with fl56 



E 



In T 



2 ln(5„) + E 



In (Tr '^' I 



(58) 



(59) 



(60) 
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Angle c 



Figure 2: F as a function of the angle c. 



and from f l50|) . we have: 



E 



2In(g„) + E 



In T. 



2c2 



b{i)i 



27r 



(61) 



For the first moment of In ( Td j , for an angle c, using Bougerol's identity 
|Bou83l IReY99] . there is the integral representation |CMY98l IVaklO) : 



E 



In (rf = 2F(c) + In (2) + cs, 



wher^: 



Fic) 



dz 



cosh (f ) 



In (sinh {cz)) 



(62) 



(63) 



and Ce ~ 0.577 denotes Euler's constant. 

For c = 27r, we have -F(27r) 3.84. In Figure [2] we plot F with respect to 
the angle c. 



^^We note that there is a simple relation between F and G: F' (^j^^ = c G(c). For 
a more complete discussion, see [VaklO| . 
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In summary, from (H7!) . (159|) . (16T|) and ( !62|) . we approximate 



-27r 



2-K 



-|21n(5„) + Q + ^i? 



l--'27r 



b(1)| 



where: 



Q = 2F(27r) + 2 In 2 + 
is a constant with F(27r) ^ 3.84, c^; ~ 0.577, and E 

Q ^ 9.54. 



b(i), 



(64) 
(65) 

0.167, thus: 
(66) 



Thus, by taking c = 2tt, for n large, the mean time E [r„] that X„(t) rotates 
around 0, is: 



o, /- N 0-08 ^ 
21n(c„) + ^ + Q 



and since: 



the MRT of Xn{t) is given by the formula: 



E [fn] ^ Y 



n 



21n^ +0.08— + Q 



(67) 



(68) 



(69) 



is 



For a long enough polymer, such that uIq » L ^ n » I, thus / = 
negligible with respect to X„(t) and from ( I3T1) X„(t) ^ X„(t). For a mean 
initial configuration Cn = E {YMt=i ^^^''^^'^) ^ obtain that the MRT of the 
polymer is given by the formula: 



E [rJ ^ 



2 In 



+ 



O.OSn 



(EiELie^''^^'^))' 



+ Q 



(70) 



Finally, using the unsealed variables (jl]), we obtain: 



E[rr, 



8D 



2 In 



0.08n 



n 



+ Q 



(71) 
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Figure 3: MRT of the free polymer end as a function of the number 
of beads n (Brownian simulations). 

Expressions (170|) and ( !7T]) show that the leading order term of the MRT 
depends on the initial configuration, however this dependence is weak. 

We consider now that the polymer is initially stretched {6k{0) = 0, V/c = (1, n)), 
hence c„ = n. Thus, from (17T1) . the MRT is approximately: 



In order to check the range of validity of formula fl72|) . we ran some 
Brownian simulations. In Figure |3l we simulated the MRT with a time step 
dt = 0.01 for n = 50 to 300 rods in steps of 5, and for each n, we took 
300 samples and averaged over all of them. The parameters we chose were 
D = 10 for the diffusion coefficient, L = 0.3 for the distance from the origin 
and Iq = 0.25 for the length of each rod, and for the initial condition we chose 
a stretched polymer located on the half line Ox, 6k{0) = 0, VA; = (l,...,n) 
(hence c„ = L + uIq), and then we computed the MRT E [r„]. In Figure HJ 
we plot both the results from Brownian simulations and the formula (172 p . 
We considered values from n = 50 to 300 rods in steps of 10, D = 10, 
L = 0.3, /o = 0.25, thus the condition ^/n Iq » L is satisfied. For each 
numerical computation, we performed 300 runs with a time step dt = 0.01. 
By comparing the numerical simulations (Fig. Hj), with the analytical formula 
for the MRT, we see an overshoot. 




ln(r2) + 0.08- + Q 



(72) 



8D 



n 
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Figure 4: Comparing the Analytical Formula and the Brownian 
simulations of the MRT. The MRT of the polymer end is depicted with 
respect to the number of beads n. 



The straight initial configuration is no restriction to the generality of our 
study: we have c„ < n and the upper bound is achieved for a straight initial 
configuration. In Figure we simulate the MRT E [r„] for both a random 
and an initially straight configuration (Brownian simulations with the same 
values for the parameters as above, after formula (172|) ). 



Uniformly distributed initial angles 



When the initial angles (^i(O), ^2(0), • • • , ^n(O)) are uniformly distributed over 
[0, 27r], by averaging over all possible initial configurations, from ([7]) and ([8]), 
we obtain: 



Cn = E 

.k=l 

hence, from ffT6l). 



E 



^(cos(e,(0)) + z sin(e,(0))) 



k=l 



0, 



E[Xr,{t)\=l + Cn = l. 



We define: 



X„(t) = — ( X„(t) -/)+! = ! + 



(73) 

(74) 
(75) 
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3.6 




Figure 5: MRT of the free polymer end as a function of the number 
of beads n for the straight and for a random initial configuration 
(Brownian simulations). 



We know that, for n large, with 0n{t) denoting the total angle of Xn{t) the 
mean time E [f„], where f„ = inf{t > 0, |<^n(^)| = 27r}, that rotates 
around 0, is: 



Finally, for a long enough polymer, such that nlo » L ^ n » I and 
Ylli=i e*^'=*-*^ >> thus I = and ^/n are negligible with respect to 
and from (175|) X„(t) ^ Xnit). Using the unsealed variables @, with Q ^ 
9.62, the MRT satisfies: 



E [r„] ^ E [f„] ^ l^g. (76) 



Remark 4.2 Formula jT/l]) or formulas F/^) and / TT^) provide the asymptotic 
expansion for the MRT when 9i G M+. In fact, 9i is a reflected Brownian 
motion in [0,27r], thus a better characterization is to estimate the MRT by 
using the probability density function for 9i in the one dimensional torus 
[0, 27r] . In the Appendix of IVakll^ . we derive this probability density function 
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n-2 

Figure 6: Brownian simulations of the MMRT as a function of the 
number of beads n — 2 for D = 10. 

and by repeating the previous calculations, we show that, e.g. formula ^7^, 
remains valid. 

4.1 The Minimum Mean First Rotation Time 

The Minimum Mean Rotation Time (MMRT) is the first time that any of 
the segments of the polymer loops around the origin, 

MMRT = min E{r„ |c„} = E [r„,„] , (77) 

where Sn is the ensemble of rods which can travel up to the origin. The 
MMRT is now a decreasing function of n. In Figure El we present some 
simulations for the MMRT as a function of n (100 simulations per time step 
dt = 0.01 with n = 4 to 15 rods, D = 10 for the diffusion constant , L = 0.3 
for the distance from the origin and Iq = 0.25 the length of each rod). 
The initial configuration is such that |<^„(0)| < 2tt — e, e.g. the straight 
initial configuration: ^^(0) = 0, V/c = (1, ...,n). In Figure[71 we present some 
Brownian simulations for the MMRT as a function of D and of L (Figure 
Wi^) and Figure [7|^b) respectively). L and Iq satisfy the rotation compatibility 
condition nlo > L. E [Tmin] decreases with D and increases with the distance 
from the origin L. It remains an open problem to compute the MMRT 
asymptotically for n large. 
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Figure 7: Brownian simulations of the MMRT as a function of the 
number of beads n — 2. (a) For several values of D with L = 0.3, (b) for 
several values of L with D = 10. 

4.2 Initial configuration in the boundary layer 

When the polymer has initially almost made a loop, we expect the MRT to 
have a different behavior. We look at this numerically, and we start with an 
initial polymer configuration in the boundary layer : 

2n-6 < \<^n{0)\ < 27r. 

where <^n is defined in eq. ([3]). The rotation of the polymer will be completed 
very fast and using Brownian simulations (100 runs per point with n = 10 
rods and D = 1 for the diffusion constant, L = 0.1 for the distance from 
the origin and Iq = 0.2 for the length of each rod), we plotted in Figure [8] 
the results showing that when the initial total angle v?n(0) tends to 27r, the 
MRT tends to zero, with the precise asymptotic remaining to be completed. 
Numerically, we postulate that there is a threshold for an initial total angle 
lv^n.(0)| = f . When |99„(0)| < |, the MRT appears to be independent from 
this angle, whereas for |y9„(0)| > ^, the MRT decreases to zero. However, 
this needs further investigation. 

5 Discussion and conclusion 

In the present paper, we studied the MRT for a planar random polymer con- 
sisting of n rods of length /q. The first end is fixed at a distance L from 
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3pi/2 



7pi/4 



2pi 



Initial Total Angle 



Figure 8: Brownian simulations of the MRT as a function of the 
initial total angle (/3n(0), where ipn is defined in ([3]). 

the origin, while the other end moves as a Brownian motion. Interestingly, 
we have shown here that the motion of the free polymer end satisfies a new 
stochastic equation (13T]) . containing a nonlinear time-dependent determinis- 
tic drift. When n is large, the limit process is an Ornstein-Uhlenbeck process, 
with different time scales for each of the two coordinates. 

We found that the MRT E [r„] actually depends on the mean initial con- 
figuration of the polymer. This result is in contrast with the one of the small 
hole th eory |WaK93l IWHK931 IHoS04l ISSHE061 ISSHOGal ISSHOGbl ISSH07| 
ISSH08] where the leading order term of the MFPT for a Brownian particle 
to reach a small hole does not depend on the initial configuration. Although 
the MRT is not falling exactly into the narrow escape problems, it is a rare 
event and the polymer completes a rotation when the free end reaches any 
point of the positive x-axis. The reason why the initial configuration survives 
in the large time asymptotic regime, is due to the dynamics of the free mov- 
ing end, approximated as a sum of i.i.d. variables, which is not Markovian. 
It leads to a process with memory. In summary, for n/o » L, the leading 
order term of the MRT is given by: 

1. for a general initial configuration: 
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where Cn = E (^^=1^^^''^^^^ , (6*^(0), 1 < A; < n) is the sequence of 
the initial angles and Q ^ 9.54, 

2. for a stretched initial configuration: 



E [Tn] ^ ^ (ln(n) + Q) 



3. for an average over uniformly distributed initial angles: 



where Q ^ 9.62. 

As we have shown, these formulas are in very good agreement with Brownian 
simulations (Fig. H]). 

After completion of this work, several questions arise naturally, namely: 

- When the polymer has already made one loop, what is the probability that 
it makes a second loop before unwrapping and in that case, what is the MRT? 

- How can we extend our study in dimension 3 and higher dimensions? See 
|H"VY11| for some first results. 



A Appendix: Proof of Theorem 13.1 



We show the convergence of the sequence {M^"'\t > 0) to a Brownian motion 
in dimension 2, with a different time scale function for each coordinate. This 
classically involves 2 steps (see e.g. |Bil68l[Bn78] or |ReY99| (Chapter XIII)): 

1. the convergence of the finite dimensional distributions, and 

2. the tightness of the sequence {M^"'\t > 0) 

a) The martingales Sj""* and Cj"-* may be written in the form of stochastic 
integrals, as: 

Ci") = -r / y2cos{Bk{s))dB,{s). 
Jo r-f 



k=l 

25 



Consequently, with < M >t denoting the quadratic variation |ReY99] of the 
martingale {Mt,t > 0), we obtain: 

(^^"^>i = -/ Vsin2(S,(s))rfs ^ / E [sm'iB^is))] ds, i7S) 
n Jq "->oo Jq 

= - / y2cos\Bk{s))ds ^ / E [cos\Bi{s))] ds, (79) 
n Jo Jo 

1 ft ft 
^^(n)^^(n)^^ = _ / Vsin(25fc(s))rfs ^ / E[smi2B,{s))]ds, 

Jo , _i n^oo Jq 



(80) 



where the classical ergodic theorem yields the a.s. convergence |Bil68] . We 

now give explicit expressions for the 3 right-hand sides. 

Consider: 

E [{exp{tB^{s))f] =E[exp{2tB,{s))] = 
from which we deduce: 

E [cos2(5i(s)) - sin2(5i(s))] = E [cos(25i(s))] = e-^\ 

and 

E[sin(25i(s))] = 0. 

Finally, we obtain: 

E[cos^(i?i(.))] =i±^, (81) 



E [sm\B,{s))\ 



l-e 



-2s 



(82) 



, . . . , iv±i 



h) The previous results allow to obtain the convergence in law for yM]:^ 
say, as n — 7- oo. Indeed, by taking f,g : R_|_ — > R, simple functions, we may 
write: 



in 



E 



^E[exp (^ S^^))], 



(83) 
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where: 



{f{u)dS^-^ + g{u)dCl^^) 



(84) 



as follows: 

E [exp {t S^^))] = E 



exp 1 z + ^ < >, 



exp 



1 



< 



.(85) 



The results obtained in part a), together with the fact that |ReY99] : 

1 



E 



exp 1 2 S^^) + - < >, 



1 



(86) 



now yield: 

E[exp {t S^^))] — >exp 



v^(^)^^+/(^) 



1 + e 



-2m 



rfu .(87) 



c) It now remains to prove the tightness |Bil78] of the distributions of the 
sequence M^^\ which follows from a classical application of Kolmogorov's 
criterion; indeed, for /3 > and C/? a positive constant, 



E 



2/3' 



< cjE\{{S^^^^)^-{S^-^)y 



+ E 



< 2cR\t-sf. 



(88) 



We refer the reader who may want more details about the arguments used to 
R,eY99) . Chapter XIII, where convergence in distribution on the canonical 



space C(M+,]R) is discussed. 
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